Precision neutron interferometric measurement of the nd coherent neutron 
scattering length and consequences for models of three-nucleon forces 
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We have performed the first high precision measurement of the coherent neutron scattering length 
of deuterium in a pure sample using neutron interferometry. We find b n d — (6.665 ± 0.004) fm 
in agreement with the world average of previous measurements using different techniques, b n d = 
(6.6730 ± 0.0045) fm. We compare the new world average for the nd coherent scattering length 
b n d = (6.669 ± 0.003) fm to calculations of the doublet and quartet scattering lengths from several 
modern nucleon-nucleon potential models with three-nucleon force (3NF) additions and show that 
almost all theories are in serious disagreement with experiment. This comparison is a more stringent 
test of the models than past comparisons with the less precisely-determined doublet scattering length 
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PACS numbers: 03.75.Dg, 21.45.+V, 25.10.+S, 25.40.Dn, 21.30.-x 
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The nature of three-nucleon forces (3NF) has been an 
especially active area of study in nuclear physics over 
the last decade and a half [3." Nonetheless, the current 
models are incomplete and have a tendency to resolve 
discrepancies in some observables at the expense of exac- 
erbating discrepancies in others 0,0]. In this Letter, we 
present a new measurement of the coherent nd scatter- 
ing length, b n d, which can be used as part of a stringent 
program of tests of NN and 3NF models. We show that 
almost all modern nucleon-nucleon (NN) potentials, in- 
cluding those with parametric 3NF adjusted to replicate 
the triton binding energy, disagree with the new world 
average value of b n( i- 

The present family of modern NN potentials — AV180], 
CD-Bonn||, and the various Nijmegen potentials^ — fit 
the extensive database of NN bound state properties and 
scattering observables to within the precision of the data. 
More recent potential models based on chiral perturba- 
tion theory predict NN properties and observables with 
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comparable accuracy 0, 0- There is growing confidence 
that the NN interaction is understood well enough that 
deviations from NN force model predictions in 3N ob- 
servables for which relativistic corrections are negligible 
can confidently be interpreted as 3NF effects. 

The simplest systems that exhibit 3NF effects are the 
bound states of 3 H and 3 He and the scattering states of 
nd and pd. Since 3 H and nd are free from long-range 
electromagnetic interactions that complicate both the- 
ory and experiment, they are the systems of choice for 
precision tests at low energy. The computational tools 
presently available to analyze these systems are believed 
to be excellent. 

The 3N interactions employed in realistic calculations 



have evolved to reconcile a growing number of discrepan- 
cies between theoretical and experimental observables 0; 
the nucleon binding energies, the nucleon-deuteron (Nd) 
doublet scattering lengths, the Nd vector analyzing 
power A y at low energies and N d sp in observables 
above the deuteron breakup threshold 11] . Although the 
trinucleon binding energy and the Nd doublet scattering 
lengths a^d have been observed to be computationally 
correlated, as shown in the well-known Phillips lines, they 
are not equivalent. The doublet scattering length has 
special importance for effective field theory (EFT) calcu- 
lations of 3N observables. A momentum cut-off parame- 
ter required in the renormalized theory can be adjusted 
to match the experimental value of 2 aAf^[l^]. which then 
forms the fundamental expansion parameter for the sys- 
tem in this channel [l3j. 

The nd bound state coherent scattering length b is the 
forward scattering amplitude due to the strong interac- 
tion from an unpolarized ensemble in the limit of zero 
neutron energy. It is simply related to the free nuclear 
doublet and quartet scattering lengths 2S+1 a n d|b| 
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where m„ and nid are the neutron and deuteron masses. 

Using the Neutron Interferometer and Optics Facility 
(NIOF) at NIST0, we have measured b n d using a pure 
T>2 gas sample. Since the quartet scattering length can be 
accurately calculated from NN potentials and is largely 
determined by deuteron properties and therefore insen- 
sitive to 3NF effects, we believe that comparison to the 
coherent nd scattering length is a more stringent test of 
NN+3NF models than comparison to the doublet scatter- 
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FIG. 1: A schematic view of the Si perfect crystal neutron 
interferometer with gas cell and quartz alignment flag. Indi- 
vidual components are discussed in the text. 

ing length. Moreover, since the effective range function in 
the quartet spin channel is a smooth function of energy, 
the quartet scattering length can be accurately extracted 
from an energy dependent phase shift analysis 0], so 
that both theoretically and experimentally knowledge 
of the bound coherent scattering length is equivalent to 
knowledge of the doublet scattering length. 

In neutron interferometry jl4j |. the phase shift of the 
neutron beam due to the optical potential of a material 
is measured. To high accuracy, the phase shift due to the 
sample is 

= (n r - l)kD eff = - WbiD eff , (2) 
i 

where n r is the real part of the index of refraction, 
(n r — 1) « 10~ 5 , k is the wavevector, and the sum is 
taken over the elemental species. To measure the bound 
coherent scattering length b to 0.05 % absolute accuracy, 
measurements of the neutron optical phase shift Ac/), the 
atom density Ni, the sample thickness D e g , the neutron 
wavelength A, and the gas purity, at the 0.02 % level are 
required. 

Measurements were performed using a perfect crystal 
silicon neutron interferometer with high phase contrast 
(80 %) and long-term phase stability (« 5° per day)0. 
The single crystal interferometer is schematically illus- 
trated in Fig. ^ A monochromatic cold neutron beam 
(E = 11.1 meV, A = 0.271 nm, A A/ A < 0.5 %) is coher- 
ently divided near point A by Bragg diffraction into two 
beams that travel along paths / and //. These beams 
are again split near points B and C and are coherently 
recombined to interfere near point D. 

The phase shift, A</>, is measured by a secondary sam- 
pling method in which a phase shifter is positioned across 



both beams. The intensities of the O and H beams as a 
function of the phase flag angle, <5, can be described by 
the following relations: 

Io(5) = A +Bcos(CaS) + AM, (3) 
I H (6) = A H + B cos (C£,(S) + A<f) + it) . 

The function £(<5) is the path length difference of the two 
beams (I and II) traversing the phase flag. The param- 
eters Ao, Ah, B, C, and Atfio are extracted from fits to 
the data. The value of A</> and its corresponding uncer- 
tainty is used to determine the phase difference between 
the two interfering beams. 

In order to minimize the effect of the phase shift due to 
the aluminum cell, the cell walls were designed to extend 
across both beam paths, producing compensating phase 
shifts. When the cell is perfectly aligned, the beams 
strike both cell compartments perpendicular to their sur- 
faces, minimizing systematics. The resulting cell phase 
shift was five times smaller than the gas phase shift, re- 
ducing the relative contribution of the cell to the total 
phase shift uncertainty to 1.7 x 10~ 4 . Alignment of the 
cell was performed by a phase shift measurement using 
a quartz plate mounted on the kinematic mount of the 
cell, similar to the procedure used in Ref. [l^ |. 

Gas was then introduced into the chamber on path J, 
while the cell on path II was evacuated. Interferograms 
with the cell first in the beam and then translated out 
of the beam were collected to determine the phase shift 
due to the gas and to the small difference in aluminum 
thickness. This phase difference is combined with mea- 
surements of N, D e ff, A, and the gas purity to extract 
b n d- 

The atom density N is determined using the ideal gas 
law with virial coefficient corrections up to the third pres- 
sure virial coefficient. For deuterium, the virial coeffi- 
cients have been measured with sufficient accuracy to de- 
termine N with a relative uncertainty of 0.001 %[lii.llfSj. 
The absolute temperature was measured using two cali- 
brated 100 fl platinum thermometers which have an ab- 
solute accuracy of 0.023 % at 300 K. The pressure was 
measured using a silicon pressure transducer capable of 
measuring the absolute pressure to better than 0.01 %. 

The gas cell was (1.0016 ± 0.0001) cm thick at an 
absolute temperature of (20.00 ± 0.05) °C. Dimensional 
changes in the cell due to the « 12 bar pressure result 
in a change in thickness at the center of the cell of less 
than 1 /zm, which amounts to a systematic effect on the 
thickness of less than 0.01 %. 

Both mass spectroscopy and Raman spectroscopy were 
employed to measure impurities in the D2 gas sample, 
with the primary contaminants being HD and D2O. The 
mole fractions of D, H, and O were measured to be 
x D = 0.99840 ± 0.00017, x H = 0.001500 ± 0.000065, and 
x = 0.000050 ± 0.000016. The expression for b nd , cor- 
rected for impurities is, b n d = [b ga s — bnXH — bo^o]/ 'xd, 
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FIG. 2: Values for the coherent scattering length on a run- 
to-run basis. 



FIG. 3: Bound coherent neutron scattering length for nd 
along with reported uncertainties. Our result for D2 value 
is consistent with the previous and the current world average. 



where b H = (-3.746 ± 0.0020) fm and b = (5.805 ± 
0.004) fm[2(j. Corrections due to the molecular state of 
the gas are not relevant within the reported accuracy of 
the measurement. 

The neutron wavelength was measured using an an- 
alyzer crystal placed in the H-beam as described in 
Ref. I2H. Rotating this crystal through both the sym- 
metric and anti-symmetric Bragg reflections allows the 
absolute Bragg angle, 6b, to be determined to high accu- 
racy, yielding A = (0.271266 ± 0.000012) nm. The stabil- 
ity of the wavelength over time was shown in a separate 
test to be better than 0.001 %. 

The value of b n d was extracted for each 42 min data set 
on a run-by-run basis (Fig. ■ These values were then 
averaged to obtain the coherent neutron scattering length 
b n d = (6.665 ± 0.004) fm. To compare this value with the 
previous world average value of b n d, we consulted the 
existing compilations of previous measurements [2fl l22l | 
and excluded all measurements which were not published 
in a referecd journal and all measurements which were 
later retracted. The remaining values were combined into 
an average with results weighted by the inverse square of 
their uncertainties in the usual manner |2fl|. The values 
are presented in Fig. The result of this evaluation is 
b nd = (6.6730 ±0.0045) fm with a reduced X 2 = 0.6. Our 
value, b n d — (6.665 ± 0.004) fm, is consistent with the 
previous world average and of comparable precision. The 
new world average value for the bound nuclear scattering 
length of deuterium is b n d — (6.669 ± 0.003) fm. 

For the forward scattering, which is the only compo- 
nent that contributes to the phase shift seen by the inter- 
ferometer, the scattering amplitudes from both the Mott- 
Schwinger and neutron-electron interactions are zero|l4|. 
The size of the scattering due to the electric polarizabil- 
ity of the neutron is less than -0.000017 fm[24|. This 
measurement is thus sensitive only to the same nuclear 
interactions that are calculated in theoretical models of 
3N scattering. 

To illustrate the impact of comparing theoretical cal- 
culations to the precision data on the coherent scattering 



TABLE I: Theoretical calculations of the nd scattering 
lengths. The rightmost column is the coherent scattering 
length. The 3NF parameters in the boldfaced potential mod- 
els have been adjusted to reproduce the triton binding energy. 



Ref. 


Potential model 


2 a nd [fm] 


4 a nd [fm] 


and [fm] 


[25] 


MT I-III 


0.70 


6.44 


4.52 


[26] 


AV14 


1.35 


6.38 


4.70 




SSCC 


1.32 


6.41 


4.71 


[27] 


RSC 


1.52 


6.302 


4.71 




RSC+TM3NF 


0.393 


6.308 


4.336 




AV14 


1.200 


6.372 


4.648 




AV14+BR3NF 


0.001 


6.378 


4.252 




RSC+TM3NF 


0.657 


6.304 


4.422 




AV14+BR3NF 


0.567 


6.380 


4.442 


[28] 


MTI-III 


0.71 


6.43 


4.52 




MTI-III 


0.702 


6.442 


4.529 




AV14 


1.196 


6.380 


4.652 


m 


AV14 


1.189 


6.379 


4.649 




AV14+TM3NF 


0.5857 


6.371 


4.443 



length, a number of modern calculations of the nd scat- 
tering lengths are shown in Table [U The dependence 
of the theoretically calculated 2 S 1 / 2 scattering length on 
the inclusion of a 3NF is clearly seen. We note that, as 
expected, none of the theories which do not incorporate 
a 3NF come close to matching the nd coherent scattering 
length. The MTI-III model, which gives fair agreement 
without a 3NF, nonetheless fails to reproduce the NN 
partial wave phase shifts. Of the potential models that 
include a 3NF, only the AV14 potential with the Brazil 
3NF[i3| and the AV14 potential with the TM 3Nf|3oT] 
are in agreement with the data as shown in Fig. 01 We 
wish to stress that the precision with which the coherent 
nd scattering length is known can, in concert with the 
triton binding energy, set tight constraints on NN po- 
tential models as well as on 3BF models. None of the 
potential models listed in Table [I] for example, include 
charge independence breaking (CIB) effects. The authors 
of a recent paper [3l| estimate that CIB effects increase 
the nd doublet scattering length by 0.19 fm. Arriving 
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at a simultaneously correct NN and 3N potential model 
requires sufficient precision in all relevant fundamental 
low-energy parameters, especially the triton binding en- 
ergy and the coherent nd scattering length. 
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FIG. 4: Theoretical calculations of the coherent scattering 
length compared with the new world average. The central 
dark band is the la confidence band and the lighter band is 
the 2a confidence band. 

We observe that none of the theories, with the excep- 
tion of the AV14 potential with a 3N force, are in agree- 
ment with the precisely-known world average coherent 
nd scattering length of a nd = (4.443 ± 0.002) fm. The 
precision of the measurement of this critical parameter is 
now such that we can distinguish between different theo- 
retical combinations of NN and 3NF models, even when 
the 3N forces have been adjusted to replicate the triton 
binding energies. Relativistic corrections to the triton 
binding energy, which were shown by the authors of [3^ 
to weaken the binding by 0.29-0.43 MeV, further reduce 
the correlation between h nd and the triton binding en- 
ergy, and provide another justification for using this pair 
of parameters to perform rigorous tests of NN and 3NF 
models. Clearly the theoretical community now has an 
incentive to systematically investigate nuclear force mod- 
els for these parameters with comparable accuracy, espe- 
cially since the precision of interferometric measurements 
of neutron-nucleus scattering can be extended by at least 
an order of magnitude. 
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